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We present our experimental investigations on the subject of dynamical nonlinearity-induced 
instabilities and of nonlinear Landau-Zener tunneling between two energy bands in a Rubidium 
Bose-Einstein condensate in an accelerated periodic potential. These two effects may be considered 
two different regimes (for small and large acceleration) of the same physical system and studied with 
the same experimental protocol. Nonlinearity introduces an asymmetry in Landau-Zener tunneling; 
as a result, tunneling from the ground state to the excited state is enhanced whereas in the opposite 
direction it is suppressed. When the acceleration is lowered, the condensate exhibits an unstable 
behaviour due to nonlinearity. We also carried out a full numerical simulation of both regimes 
integrating the full Gross-Pitaevskii equation; for the Landau-Zener effect we also used a simple 
two-level model. In both cases we found good agreement with the experimental results. 

PACS numbers: PACS numbor(s): 03.65.Xp, 03.75.Lm 
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I. INTRODUCTION 



Cold atoms and, more recently, Bose-Einstein conden- 
sates (BECs) in optical lattices have attracted increas- 
ing interest since their first realization In particu- 
lar, the formal similarity between the wavefunction of 
a BEC inside the periodic potential of an optical lat- 
tice and electrons in a crystal lattice have triggered the- 
oretical and experimental efforts alike. Many phenom- 
ena from condensed matter physics, such as Bloch os- 
cillations and Landau-Zener tunneling have since been 
shown to be observable also in optical lattices 0, 0, Q . 
In a recent experiment, a BEC in an optical lattice even 
made possible the observation of a quantum phase tran- 
sition that had, up to then, only been theoretically pre- 
dicted for condensed matter systems 0. However, an 
important difference between electrons in a crystal lat- 
tice and a BEC inside the periodic potential of an op- 
tical lattice is the strength of the self interaction and 
hence the magnitude of the nonlinearity of the system. 
Electrons in a metal are almost noninteracting whereas 
atoms inside a BEC interact strongly. A perturbation 
approach is appropriate in the former case while in the 
latter the full nonlinearity must be taken into account. 
From this feature new physics is expected. Most exper- 
iments to date have been carried out in the regime of 
shallow lattice depth, for which the system is well de- 
scribed by the mean field Gross-Pitaevskii equation with 
a periodic potential. Moreover, the nonlinearity induced 
by the mean-field of the condensate has been shown, both 
theoretically and experimentally, to give rise to instabili- 
ties 6, 7, 8, 9, jr., J^2, 13, 14J in certain regions of the 
Brillouin zone. These instabilities are not present in the 
corresponding linear system, i.e. the electron system. 



In this paper we review and summarize our experi- 
mental and theoretical results on the subject of nonlinear 
Landau-Zener tunneling and nonlinearity-induced insta- 
bilities in a Bose-Einstein condensate interacting with an 
external periodic potential. These two phenomena rep- 
resent the most dramatic manifestations of nonlinearity 
in two different regimes of the system. In order to study 
these phenomena we have used a single experimental pro- 
cedure. The underlying idea is to linearly scan the Bril- 
louin zone, by applying a constant acceleration to the 
periodic potential, and cross the band edge. Then the 
condensate is released and an absorption picture of the 
condensate is taken after a time of flight, reflecting the 
momentum distribution at the time of release. By vary- 
ing the nonlinearity of the system with a flxed (large) 
acceleration, we can study the nonlinear contributions to 
the Landau-Zener effect. We shall denote this regime as 
the "Landau-Zener" regime. On the contrary, by varying 
the acceleration from very small to intermediate values 
with a fixed nonlinearity, we can study the stability of 
the condensate and the effects of nonlinearity on the dy- 
namics. We will refer to this regime as the "instability" 
regime. 

This paper is organized as follows. After describing 
our theoretical approach in section^ we explain our ex- 
perimental techniques in section UllI Section Hvl presents 
a discussion of our results on the Landau-Zener tunnel- 
ing, and the experimental and conceptual difficulties en- 
countered in obtaining them. For the interpretation of 
the nonlinear Landau-Zener tunneling we re-examine and 
critically compare the effective potential concept, intro- 
duced into previous investigations, with our present re- 
sults. Furthermore we interpret the asymmetry in the 
nonlinear Landau-Zener effect on the basis of different 
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chemical potentials calculated for the ground band and 
for the first excited band. Section[V]discusses our experi- 
mental and theoretical results on the condensate instabil- 
ities. Finally our conclusions and perspectives for future 
developments and improvements are given in section IVTl 



II. THEORY 

The motion of a Bose-Einstein condensate in an ac- 
celerated ID optical lattice is described by the Gross- 
Pitaevskii equation 
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FIG. 1: Schematic representation of the band structure for 
the chemical potential /i versus the quasimomentum q in an 
optical lattice (s = 2) and LZ-tunneling (ground to excited 
band (a) and excited to ground band (b)). When the BEG 
is accelerated across the edge of the Brillouin zone (BZ) at 
quasimomentum qs = hkh, LZ tunneling can occur. Further 
acceleration will result in the condensate part in the upper 
level, to undergo LZ tunneling to higher bands with a large 
probability (due to the smaller gaps between higher bands) 
being that part essentially unaffected by the lattice. After the 
first crossing of the BZ edge, we increase the lattice depth and 
decrease the lattice acceleration thus reducing the tunneling 
rate from the ground band to the excited band to much lower 
values at successive BZ edge crossings. 
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where M is the atomic mass, ~ ir/d is the optical 
lattice wavenumber with optical lattice step d. Vq is the 
periodic potential depth, Eu = fi?k\/2M is the recoil 
energy. We introduce the dimensionless parameter s 



Va = sEr 



(2) 



denoting the lattice depth in units of the recoil energy. 
The s-wave scattering length determines the nonlin- 
earity of the system, with the two-body coupling constant 
given by Airti^as/M. Equation is written in the co- 
moving frame of the lattice, so the inertial force Maj^ ap- 
pears as a momentum modification. The wavefunction tp 
is normalized to the total number of atoms in the conden- 
sate and we define no as the average uniform atomic den- 
sity. By defining the dimensionless quantities x = 2kLX, 
t = 8ERt/H, ij = tp/y/T^;, V = s/16, a = MaL/WEakL, 
qb = 1/2. The nonlincarity is characterized through the 
parameter [l5l | 



(3) 



Therefore eq. ^ is cast in the following form p^ . lrRIT^ : 



dx 



at] ii + vcos{x)ij + C\ip\^ (4) 



where we have replaced x with x, etc. In the neigh- 
borhood of the Brillouin zone edge, at quasimomentum 
q = qb , approximate the wave function by a superpo- 
sition of two plane waves with complex coefficients (the 
two level model illustrated in ^^), assuming that only 
the ground state and the first excited state are popu- 
lated 0. We then substitute in eq. I@J 



(5) 



with |a(t)p + [6(t)P — 1. Comparing the coefficients of 
^iqx g^j^j ^i{q-i)x ^ linearizing the kinetic terms and drop- 
ping the irrelevant constant energy 1/8 + C[l + (|ap + 
|6p)/2], eq. Q assumes the form 
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where (t^ (i = 1,2,3) are the Pauh matrices. Each solu- 
tion of eq. © has an associated conserved energy e[i/'] 



e[-0] ^ at < V'l ylV' > +v < V'l yl''/' > + 



C(< ^lylV- >)^ with IV >E 



(7) 



It must be stressed that the energy and the hamilto- 
nian eigenvalue do not coincide, since we are dealing with 
a nonlinear system. The energy is the quantity that is 
conserved along the system trajectories while the hamil- 
tonian eigenvalue is not. In nonlinear systems the hamil- 
tonian eigenvalue is usually called the chemical potential 
(indicated with /i) . The connection between energy e and 
chemical potential /z for a condensate within an optical 
lattice was discussed in ^ . Those quantities are equiv- 
alent in the case of negligible atomic interactions, i.e. 
C ~ 0. The Bloch bands may describe either the energy 
or the chemical potential as a quasimomentum function. 
In the following we choose to consider the Bloch bands 
as the curves representing the chemical potential /i as a 
function of the quasimomentum (fig.0. 

The adiabatic Bloch bands of eq. © have a swallow 
tail structure (i.e. they develop a loop and become multi- 
valued) at the edge of the Brillouin zone for C > w 0, 
1^122, 113, a regime not explored in our experiments. 

In fig. EI a) we plot the band gap at the BZ edge be- 
tween the lower band and the first excited band in the 
nonlinear case, calculated in regimes that are relevant 
for the Landau-Zener experiments. We defined the gap 
in the non linear case as the chemical potential difference 
between the two bands at the edge of the BZ. However 
a definition of this kind is manifold: we define A/x as 
the difference in chemical potential between the first ex- 
cited and ground band, both calculated self-consistently 
assuming the total density of the condensate in each of 
them; moreover we define A/j,io as the chemical poten- 
tial difference between the first excited and ground band, 
both consistent with the total density in the ground band 
(i.e. using the ground state wavefunction in the mean 
field term) and A/ioi as the chemical potential difference 
between the first excited and ground band, both consis- 
tent with the total density in the first excited band (i.e. 
using the excited state wavefunction in the mean field 
term) . In order to describe a transition from the ground 
band to the upper one, we believe that the gap A/iio is 
better suited, and conversely the gap A/xgi should better 
describe the tunneling in the opposite direction. 



III. EXPERIMENTAL SETUP 

Our experimental apparatus for creating BECs of ^^Rb 
atoms was described in |24j . The main feature of our 
apparatus relevant for the present work is the triaxial 
time-averaged orbiting potential (TOP) trap with trap- 
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FIG. 2: In a), the energy gaps between the ground and the 
excited Bloch bands calculated at the BZ edge for different 
values of the nonlinear parameter C are shown. The contin- 
uous line represents the energy gap A/i while the dashed and 
dot-dashed lines shows the energy gaps A/^io and A/^oi as 
explained in the text. The dotted lines are eq. 11211 for both 
tunneling directions. In b), the Landau-Zener tunneling rate 
r calculated from the energy gaps in a) and from eq. 11211 . with 
the same conventions for the linestyles as in a) , is shown. The 
calculations were done for a lattice depth s = 2.2. 



trap is, therefore, almost isotropic. The optical lattice 
is created by two laser beams with parallel linear polar- 
izations and wavelength A, as described in 25]. The two 
beams are derived from the first diffraction orders of two 
acousto-optic modulators that are phase-locked but with 
independent frequencies, allowing us to introduce a fre- 
quency difference /S.v between them. The resulting peri- 
odic potential has a variable lattice constant d depending 
on the intersecting angle between the two laser beams. 
The smallest lattice constant, 0.39^m, is obtained in the 
counterpropagating configuration and can be increased 
up to 1.2/xm when the two lasers intersect at about 38 
degrees. The depth Vq of the periodic potential (depend- 
ing on the laser intensity and detuning from the atomic 
resonance of the rubidium atoms) can be varied from 
Eji up to approximatively 3 Eji. It must be noted that 
Eji is the true recoil energy and depends on the angle at 
which the two laser beams intersect. The beams are de- 
tuned to the red side of the rubidium atomic resonance 
by 30 GHz. In this way, a periodic potential with lat- 
tice recoil energy En/h = 455 Hz is created. In addition, 
by linearly chirping the frequency difference Av, the lat- 
tice is accelerated with = rf^^. In our experiments. 
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we used accelerations ranging from = 0.3 ms~^ to 
= 5ms~^. 

The experimental protocol for 'moving' the condensate 
across the Brillouin zone is as follows. After creating 
BECs with roughly 10^ atoms, we adiabatically relax the 
magnetic trap frequency to Vx — 42 Hz. Thereafter, the 
intensity of the lattice beams is ramped up from Er to 
a value corresponding to a lattice depth of approxima- 
tively 2 En. Once the final lattice depth is reached, the 
lattice is accelerated for a time t. Finally, both the mag- 
netic trap and the optical lattice are switched off, and 
the condensate is observed by absorption imaging after a 
time-of- flight of 21 ms. 



IV. LANDAU-ZENER TUNNELING 
A. Modelling 

Evaluating the transition probability in the adiabatic 
approximation for the transition from an initial state to 
a final one separated by an energy gap, we find the linear 
LZ formula for the tunneling probability r 



— p 2q 



(8) 



expressing the state population changes in terms of the 
rate a at which the diagonal terms of the linear Hamil- 
tonian change their value, and of the off-diagonal inter- 
action strength v It must be noted that since we 
are considering two states only, v represents the energy 
gap, too. The transition probability is symmetric in the 
linear case, i.e. the tunneling rate from the lower level to 
the upper one is the same as in the opposite direction. 

In the nonlinear regime, as the nonlinearity parameter 
C grows, the lower to upper tunneling probability grows 
as well until an adiabaticity breakdown occurs aX C — v 
where the swallow tail structure appears 18] . The upper 
to lower tunneling probability, on the other hand, de- 
creases with increasing nonlinearity |27j . We derive the 
tunneling rate from the numerical integration of eq. 
In fig. 13 we plot the lower to upper tunneling rate (ini- 
tial (a, 6) — (1,0) in eq. JSJ) and the upper to lower 
tunneling rate (initial (a, 6) = (0,1) in eq. jsj) of the 
Bose-Einstein condensate as a function of the nonlinear 
parameter C. We see that for C — the rate is the same 
for both tunneling directions whereas for C 7^ the two 
rates are different. We confirm the presence of tunneling 
asymmetry by integrating eq. ^ directly (taking into ac- 
count the full experimental protocol), finding qualitative 
agreement with the prediction of the two-state model. 

From an analytical point of view the nonlinear regime 
is interpreted straightforwardly by writing eq. ^ as 
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FIG. 3: LZ tunneling rate r as a function of the nonlinear 
parameter C. The dashed and continuous lines represent the 
results of the two-level model of eq. ||SJ for a transition from 
the lower level to the upper one and vice versa, respectively. 
The open and filled symbols are the results of the numeri- 
cal integration of eq. Q taking into account the presence of 
the harmonic trap and simulating the full experimental proce- 
dure described in section UTTI Transition rates were evaluated 
by means of the Fourier transform method. Results are for 
V — 0.1375 corresponding to s = 2.2 and a — 0.03636 corre- 
sponding to ul ~ 2.925m s~^. C was varied in (a) by varying 
the trap frequency, and in (b) by varying the atomic scatter- 
ing length Us, the harmonic trap frequency being fixed at 20 
Hz. 



This different form leads to the same condensate dynam- 
ics since the difference between the matrices in eq. Q and 
eq. © times the (a, b) vector is equivalent to the iden- 
tity operator. The difference is then an offset in the en- 
ergy of the system. In eq. © we identify an off-diagonal 
term \v + 2C a*b\ acting as an effective potential Ve/f; the 
modulus is needed since a*b is complex. The off-diagonal 
scalar product between the two states a*b may be eval- 
uated using the adiabatic approximation technique [2^ 



1*6= -- 



2VaH^ + 



(10) 



for a transition from the lower state to the upper state. 
For tunneling in the opposite direction, a*b simply 
changes sign. The explicit expression for v^ff is 
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Within the spirit of the adiabatic approximation we put 
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a = obtaining 



^eff =V\|l± — + -^^V^l±- 



{12) 



where the upper and lower signs corresponds to initial 
conditions of excited/ground states. For the ground en- 
ergy band the following effective potential was intro- 
duced by Choi and Niu and experimentally tested 
by Morsch et al. 



H-4C 



(13) 



For small C values we modify the LZ formula of eq. ^ 
to include nonlinear corrections, replacing the potential 
V by the effective potential Vef / 0| • 



□ 


■ PRL 2003 


o 


• PRL 2001 



0.8 




0.00 0.05 0.10 0.15 0.20 0.25 



FIG. 4: Ratio between the effective potential He// and the 
bare potential w as a function of the nonlinear parameter C. 
The dotted lines are eq. (I12II for both tunneling directions. 
The solid line is eq. 1131 for tunneling from ground to excited 
band. The shaded regions are presented in the text. The 
squared marks are the experimental values calculated using 
the LZ tunneling rate of eq. Q and data from 16] for a tran- 
sition from the lower band to the upper band (open marks) 
and in the opposite direction (filled marks). The circular 
marks are data from Q with optical lattice step d — 1.18 /^m 
(filled marks) and d = 0.39 /xm (open marks), as explained 
in section UTTl All calculations and experimental data are for 
V — 0.1375 corresponding to s = 2.2 and a = 0.03636 corre- 
sponding to a_L ~ 2.925ms~^. 



Both expressions H12(l and assume to use an effec- 
tive potential taking into account the nonlinearity of the 
system. This assumption has some validity limits for the 
potential v and/or the nonlinearity strength C. Equa- 
tion H13(l . derived using perturbation theory, is expected 
to hold as long as the condensate density is nearly uni- 
form (i.e. Veff ^ 1). This condition may be realized with 
either a weak external potential or a strong atomic inter- 
action 01 . Because we consider a condensate weakly 
perturbed by the presence of the periodic potential, we 
expect eq. (|13|) to hold in a neighborhood of the center of 



the band where the effects of the periodic potential are 
weakest and the band resembles the free particle disper- 
sion curve. On the other hand we developed the two level 
model by approximating the wave function at the band 
edge; hence we expect eq. (|12|) to hold in a neighborhood 
of the band edge. 

In order to test our intuitions we performed a numer- 
ical variational solution of eq. 10} based on the Math- 
ieu functions, which are the exact solutions of the linear 
problem (i.e. eq. with C = and substituting the 
time derivative with the chemical potential). We evalu- 
ated the chemical potential using the Mathieu functions 
and we minimized this value by varying the parameter in 
the Mathieu function corresponding to the lattice poten- 
tial strength, thus obtaining an effective potential. In the 
effective potential description, this is the best approxima- 
tion of the exact solution of the nonlinear problem. We 
repeated this calculation for different values of the atomic 
quasimomentum q, from q — to the band edge q — 1/2. 
Our findings are summarized in fig.^jwhere the lower and 
upper shaded areas correspond to intermediate values 
of the quasimomentum in both the ground band (lower 
shaded area) and excited band (upper shaded area). The 
effective potential at g = corresponds to the upper edge 
of the lower shaded area (ground band) and to the lower 
edge of the upper shaded area (excited band). Figure^ 
also reports the results of eq. ifT^ and eq. ((T^ . Equa- 
tion (|13() is indistinguishable from the upper edge of the 
lower shaded area, hence we infer that it can only hold 
for the lower band exactly at quasimomentum q — 0, 
in agreement with our previous analysis. On the other 
hand, because eq. H12() reproduces very well the lower 
edge of the lower shaded area and the upper edge of the 
upper shaded area, we confirm its validity range to be 
limited to a neighborhood of the band edge. 

Equations H12|l and H13() define a different role of the 
nonlinearity. Equation H13|l states that the ratio Veff/v 
does not depend on the potential magnitude v, but only 
on the magnitude of the nonlinear parameter C, whereas 
eq. H12I) predicts a deviation from the linear case accord- 
ing to the ratio between the magnitude of the nonlin- 
earity and the bare potential strength. In the present 
description, the energy scale determined by the bare po- 
tential strength v appears to be irrelevant at the center 
of the band; the only relevant energy scale is fixed by the 
mean field interaction strength C. On the contrary, the 
energy scale determined by v acquires more and more sig- 
nificance toward the band edge. At the band edge the two 
energy scales have the same relevance and Vef / /v only de- 
pends on the ratio of them. The swallow tail threshold 
C/v — 1 represents a critical value for the system since 
the effective potential (|12ll vanishes for a transition from 
the lower band to the upper band while it doubles for 
a transitions in the opposite direction. Therefore when 
C/u > 1, it is no longer possible to interpret the nonlin- 
earity through an effective potential. 

The difference in the tunneling rates in the two direc- 
tions may be also derived from the difference between the 
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energy gaps A/iio and A/xqi (fig-EJ^) using eq. JH)) to link 
the transition rates to the energy gaps (fig-EIb)). As a 
consequence, the tunnehng rate is enhanced in one case 
and it is suppressed in the other one. These tunnehng 
rates are in very good agreement with those predicted by 
the effective potential of eq. p2|l . 



B. Experimental observations 

Landau-Zener tunneling between the two lowest en- 
ergy bands of a condensate inside an optical lattice was 
investigated in the following way (see fig. Initially, 
the condensate was loaded adiabatically into one of the 
two bands, either in the ground band or into the ex- 
cited band. Subsequently, the lattice was accelerated in 
such a way that, at the BZ edge, a finite probability for 
tunneling into the other band resulted. After the tun- 
neling event, the two bands had populations reflecting 
the Landau-Zener tunneling rate. In order to experimen- 
tally determine the number of atoms in the two bands, 
we then increased the lattice depth and decreased the ac- 
celeration. In this way, successive crossings of the band 
edge resulted in a much reduced Landau-Zener tunnel- 
ing probability between the ground band and the first 
excited band (of order a few percent), as illustrated in 
fig-D The fraction of the condensate that populated the 
ground band after the first tunneling event, therefore, re- 
mained in that band, whereas the population of the first 
excited band underwent tunneling to the second excited 
band with a large probability (around 90 percent) as the 
gap between these two bands is smaller than the gap be- 
tween the two lowest bands. Once the atoms underwent 
tunneling into the second excited band, they essentially 
behaved as free particles. For both tunneling directions, 
the tunneling rate r was derived from the ratio between 
the number of atoms experiencing the tunneling process 
and the total number of atoms measured from the ab- 
sorption picture. 

In order to verify the experimental procedure for mea- 
suring the tuneling rate in both directions and also to 
verify that in the linear regime the Landau-Zener tun- 
neling was symmetric, we measured the two tunneling 
rates as a function of the lattice depth for a conden- 
sate in a weak magnetic trap and hence a small value of 
the interaction parameter C. In this case, both tunnel- 
ing rates were essentially the same and agreed well with 
the linear Landau-Zener prediction, as reported in |16| . 
By contrast, when C was increased, the two tunneling 
rates began to differ. For instance for the parameter 
C = 0.095(9), the measured tunneling rate from the 
ground to excited state was r = 0.72 ± 0.10, whereas 
we measured r = 0.37 ± 0.05 in the opposite direction, 
proving the tunneling asymmetry. 

The effective potential description allows us to present 
within a unified picture the nonlinear Landau-Zener tun- 
neling rates measur ed in refs. as plotted in fig.H 
together with the theoretical predictions. We derive a 



qualitative agreement with the theoretical predictions of 
the non-linear Landau-Zener model, whereas quantita- 
tively there are significant deviations. We believe these 
to be partly due to experimental imperfections. In partic- 
ular, the sloshing (dipolar oscillations) of the condensate 
inside the magnetic trap can lead to the condensate not 
being prepared purely in one band due to non-adiabatic 
mixing of the bands if the initial quasimomentum is too 
close to a band-gap. Furthermore the amplitude of the 
shaded areas in fig. ^ points out a strong dependence 
of the effective potential on the quasimomentum. If the 
condensate trapped in the periodic potential has a fi- 
nite extension in momentum space, the overall effect of 
the effective potential on the tunneling probability is not 
restricted to a single quasimomentum but represents a 
mean effect over all the quasimomenta of the condensate. 
These effects could be responsible for the experimental 
points not falling exactly on the lines corresponding to 
the effective potential evaluated at the band edge. In any 
case the experimental data must fall within the shaded 
areas corresponding to the band in which the condensate 
was loaded, as in fig. ^ 

Because of the elastic force of the magnetic harmonic 
trap, it is important in the experiment not to drag the 
condensate too far from the rest position. If the conden- 
sate is dragged too far, the dragged part starts to feel the 
restoring force due to the harmonic potential and hence 
does not feel a constant force anymore. For this rea- 
son it is not possible to study large C values by varying 
only the harmonic trap frequency. For our experimental 
parameters, C « 0.11 was found to be the largest accept- 
able value, corresponding to a harmonic trap frequency 
of about 50 Hz. Furthermore, a numerical simulation of 
the experiment showed that for large values of C, for 
which the magnetic trap frequency was large, the mea- 
sured tunneling rates were significantly modified by the 
presence of the trap. However, we verified in the simu- 
lation that when C was varied without varying the trap 
frequency, the asymmetric tunneling effect persisted. In 
future experiments, one might study large C values by 
increasing the atomic density in the condensate by using 
an additional optical trap, in order to increase the radial 
trapping frequency or, alternatively, by using Feshbach 
resonances to vary the atomic scattering legth ag. The 
results of a numerical simulation using the latter method 
are reported in fig. 13 (b). 



V. INSTABILITIES 

A. Experimental observations 

In the previous experiment we investigated the depen- 
dence of the tunneling probability as a function of the 
nonlinearity magnitude with fixed acceleration. We have 
also investigated the stability of the condensate as a func- 
tion of the acceleration with (roughly) fixed nonlinear- 
ity [l2l |. When the condensate acquired a quasimomen- 
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turn close to the band edge, the unstable solutions of 
eq. ^ grow exponentially in time, leading to a loss of 
phase coherence of the condensate along the direction 
of the optical lattice. In our experiment, the time the 
condensate spent in the 'critical region' where unstable 
solutions existed, was varied through the lattice acceler- 
ation. When the acceleration was small the condensate 
moved across the Brillouin zone more slowly and hence 
the growth of the unstable modes became more im- 
portant. Figs.Et and Eli show typical integrated profiles 
of the interference pattern obtained for a lattice accel- 
eration ol — 0.3 ms~^. Here, the condensate reached 
the same point close to the Brillouin zone edge as in 
Figs. audita, but because of the longer time it spent 
in the unstable region, the interference pattern was al- 
most completely washed out. It is also evident that the 
radial expansion of the condensate was considerably en- 
hanced when the Brillouin zone was scanned with a small 
acceleration. 
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FIG. 5: Integrated longitudinal and transverse profiles of the 
interference pattern of a condensate released from an optical 
lattice after acceleration to a quasimomentum q « O.Ogs and 
a subsequent time-of-flight of 21 ms. The Gaussian width wr 
of the transverse width is marked in (b) and (d). In (a) and 
(b), the acceleration ul was 5ms~^, whereas in (c) and (d) 
aL = 0.3 ms~^. In (a) and (c), the horizontal axis has been 
rescaled in units of recoil momenta. Note the different vertical 
axis scales (by a factor 4) for the upper and lower graphs. The 
total number of atoms was measured to be the same in both 



In order to characterize more quantitatively our ex- 
perimental findings on the instability, we defined two ob- 
servables for the timc-of-fiight interference pattern. By 
integrating the profile in a direction perpendicular to the 
optical lattice direction, we obtained a two-peaked curve 
(see fig. El (a)) for which we defined a visibility V (in anal- 
ogy to spectroscopy) reflecting the phase coherence of the 
condensate. V is close to 1 for perfect coherence, whereas 
V — > for an incoherent condensate. The second ob- 
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FIG. 6: Visibility V and radial width wr as a function of the 
final quasimomentum (in units of gs) for different accelera- 
tions. As the acceleration was lowered, instabilities close to 
the BZ edge decreased the visibility and increased the radial 
width. For comparison, in each graph the (linear) fits to V 
and Wr for the a_L = 5ms~^ data are included. 



servable we defined, was the width wr of a Gaussian fit 
to the interference pattern integrated along the lattice 
direction over the extent of one of the peaks (see fig. 
(b) and (d)). 

In W6 measured V and wr as a function of the final 
quasimomentum value reached for different values of the 
acceleration. FigurejBlshows clearly that for large acceler- 
ations, both V and wr remained reasonably stable when 
the edge of the Brillouin zone is crossed. In contrast, for 
ttL — 0.3 ms~^ one sees a drastic change in both quanti- 
ties as the quasimomentum approaches the value qs ■ For 
those accelerations, the condensate spent a sufficiently 
long time in the BZ unstable region and hence lost its 
phase coherence, resulting in a sharp drop of the visibil- 
ity. At the same time, the radial width of the interference 
pattern increased. This increase is evidence for an insta- 
bility in the transverse directions. For ql = 0.3ms~^, 
the interference patterns for quasimomenta larger than 
unity were so diffuse that it was not possible to measure 
either the visibility or the radial width in a meaningful 
way. 

In [l^ we estimated an instability growth rate of 
lO'^s^^ from the time spent by the condensate in the 
unstable region of the Brillouin zone. This experimental 
value agrees with the theoretical prediction of 2000 s^^ 
estimated by Niu 0. Fallani et al. 13] carried out an 
experiment for operating conditions similar to those of 
our investigations, but measured the effect of the insta- 
bility through the atom loss rate from the condensate. 
They reported atom loss rates in the range 40 — 60 s~^. 
Although it is not clear how to link the atom loss rate 
to the instability growth rate, the qualitative agreement 
(i.e. the dependence of both quantities on the quasimo- 
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mentum) between their experimental results and their 
theory is good. Their experiment was subsequently sim- 
ulated by Modugno et al. ^ who predicted growth rates 
of the unstable modes in the range 2000 — 3000 s~^. In 
that work the growth rate is expressed in units of the 
transverse magnetic trap frequency. This choice could 
be interpreted as meaning that the instability depends 
on the radial modes. It has, however, been clarified by 
one of the authors of that paper [sOl that this is not the 
case and that the choice of the transverse frequency as 
the frequency unit was made purely for reasons of conve- 
nience. In fact, the dependence of the growth rate on the 
transverse frequency needs to be studied in more depth. 
We should point out that all the theoretical analyses pre- 
dict the occurrence of the instabilities at the microscopic 
level only, while the experimental observations reflect the 
macroscopic changes of the condensate features, for in- 
stance the visibility in our case or the atom loss rate in 
the case of 111. 




T r 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

q/qe 

FIG. 7: Results of a one-dimensional numerical simulation for 
the visibility V versus the condensate final quasimomentum 
q in the conditions of the experiment for acceleration a — 
0.3 ms~^ and for different values of the nonlinear parameter 
C. The open squares, circles and triangles correspond to C = 
0.008 (the value for our experiment), C = 0.004 and C = 0, 
respectively. The dashed lines connect the theoretical points 
to guide the eye. The closed symbols are the experimental 
values of the visibility as reported in fig.|S|for a — 0.3ms~^. 



B. Modelling 

We compared our experimental results to a simple 1- 
D numerical simulation. Figure [7| shows the results of 
a numerical integration of the one-dimensional Gross- 
Pitaevskii equation with the parameters of our experi- 
ment. The visibility was calculated in the same way as 
was done for the experimental interference patterns. It is 
clear from this simulation that it is, indeed, the nonlin- 
earity that is responsible for the instability at the edge 
of the Brillouin zone. When C was set to in the nu- 
merical simulation, the visibility remains unaltered when 
the BZ edge is crossed, whereas for finite values of C the 



visibility decreases as the quasimomentum qb was ap- 
proached. Furthermore, the larger the value of C, the 
more pronounced was the decrease in visibility near the 
band edge. For C = 0.008, corresponding to the value re- 
alized in our experiment, the onset of the instability was 
located just below a quasimomentum of 0.8qB- Experi- 
mentally, we found that the visibility started decreasing 
consistently beyond a quasimomentum of w 0.6 — 0.7 qb, 
agreeing reasonably well with the results of the simula- 
tion. The presence of experimental points with visibil- 
ity less than unity in the quasimomentum region below 
O.GqB can be explained by considering the initial sloshing 
of the condensate inside the harmonic trap. This oscilla- 
tion introduced a sensitive error in the quasimomentum 
determination which was estimated to be of the order of 
0.3gs. 

The determination of the quasimomentum correspond- 
ing to the onset of instability is a subject that has been 
examined in the literature, and in a recent experimen- 
tal investigation the quasimomentum scanning across the 
Brillouin zone was stopped at different values in order to 
verify the instability growth jl^. Our measured values 
for the instability onset are in agreement with those mea- 
sured in that reference. 



VI. CONCLUSIONS 

We have numerically simulated and experimentally 
studied the dynamics of a Bose Einstein condensate in- 
side a periodic potential in two different regimes. In the 
Landau-Zener regime, we investigated the tunneling be- 
tween two energy bands in a periodic potential and found 
that, in the presence of a nonlinear interaction term, an 
asymmetry in the tunneling rates arises. Experimentally, 
we measured these tunneling rates for different values of 
the interaction parameter and found qualitative agree- 
ment with the simulations. In the instability regime, 
we studied the stability of a BEC in the vicinity of the 
band edge, finding good agreement between experimen- 
tal results and the theoretical expectation of unstable be- 
havior. These observations confirmed that Bose-Einstein 
condensates may be used to simulate a variety of non- 
linear physics configurations. Future experiments could 
probe the complicated and time-dependent tunneling be- 
haviour due to the changing tunneling rate for multiple 
crossings of the zone edge. 

To conclude, we note that the phenomenon of asym- 
metric tunneling should be a rather general feature of 
quantum systems exhibiting a nonlinearity. For instance, 
calculating the energy shift due to a nonlinearity for two 
adjacent levels of a harmonic oscillator, one finds that 
both levels are shifted upwards in energy, the shift being 
proportional to the population of the respective level. 
The energy difference between the levels, therefore, de- 
creases if only the lower state is populated and increases if 
all the population is in the upper level. Furthermore the 
asymmetric Landau-Zener tunneling rate was applied to 
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intepret the photoassociation of a Bose-Einstein conden- tional to the rate at which the resonance is crossed, 
sate , with the surprising result that at small cross- We acknowledge financial contributions from the "Pro- 
ing rates the no-transition probability is directly proper- getto MIUR COFIN-2004" . 
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